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We discuss localized ground states of Bose-Einstein condensates in optical lattices with attractive 
and repulsive three-body interactions in the framework of a quintic nonlinear Schrodinger equation 
which extends the Gross-Pitaevskii equation to the one dimensional case. We use both a variational 
method and a self-consistent approach to show the existence of unstable localized excitations which 
are similar to Townes solitons of the cubic nonlinear Schrodinger equation in two dimensions. These 
solutions are shown to be located in the forbidden zones of the band structure, very close to the 
band edges, separating decaying states from stable localized ones (gap-solitons) fully characterizing 
their delocalizing transition. In this context usual gap solitons appear as a mechanism for arresting 
collapse in low dimensional BEC in optical lattices with attractive real three-body interaction. The 
influence of the imaginary part of the three-body interaction, leading to dissipative effects on gap 
solitons and the effect of atoms feeding from the thermal cloud are also discussed. These results 
may be of interest for both BEC in atomic chip and Tonks-Girardeau gas in optical lattices. 
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I. INTRODUCTION 



Bose-Einstein condensates (BEC) in periodic poten- 
tials are presently receiving a great deal of theoretical 
and experimental interest due to the possibility to ex- 
plore a whole class of phenomena ranging from Bloch os- 
cillations Q, , Landau-Zener tunneling Q and solitons 
to quantum phase transitions of the Mott insulator 
type [5j . In the mean field approximation these systems 
are described by the Gross-Pitaevskii equation (GPE) 
which is a cubic nonlinear Schrodinger (NLS) equation 
with periodic potential in which the cubic nonlinearity 
models the two-body interatomic interactions appropri- 
ate for dilute gases. The presence of the optical lattice 
(OL) allows interesting localization phenomena such as 
the formation of gap-solitons, i.e. localized states with 
energies comprised in the gaps of the band structure of 
the underlying linear periodic problem, both for attrac- 
tive and repulsive interactions [a, 1?J ■ This is quite re- 
markable because it is known that for two (2D) and three 
dimensional (3D) NLS, in absence of periodic potential, 
soliton solutions do not exist and for attractive interac- 
tions the phenomenon of collapse in finite time appears. 
In this last case only one localized solution is possible, 
the so called Townes soliton, existing for a single value 
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of the norm (number of atoms). This solution however 
is unstable against norm variations i.e. it collapses in a 
finite time for norms above the critical value and decays 
into the uniform background for norms below the critical 
value. The presence of the OL allows to expand the range 
of existence of localized solutions from a single value to 
a whole interval of norms below which the delocalizing 
transition occurs and above which collapse takes place. 

A description based on the GPE with cubic nonlinear- 
ity, however, is adequate only at low densities. For higher 
densities the three-body interactions start to play a role 
and a description based on two-body interactions is no 
longer sufficient. H |j| El El- To this regard we recall 
that condensate densities limited by three-body inelastic 
collisions have indeed been achieved. In these regime 
a more accurate treatment of the mean-field energetics of 
a dense condensate will need to account for both two- and 
three-body elastic collisions. 

The contribution of the three-body interactions can be 
enhanced by detuning to zero the cubic two body term by 
means of Feschbach resonances , this leading to a peri- 
odic NLS equation with quintic nonlinearity. For appro- 
priate values of density and scattering lengths, however, 
the three-body collisions could largely dominate the two- 
particle contributions even in a very dilute regime. This 
occurs when the so called Efimov effect becomes pos- 
sible i.e. when the two-body scattering length becomes 
much larger than the effective two body interaction ra- 
dius (this usually occurring near a two-body resonance). 
In this case a very large number of three-body bound 
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states can be formed in the system and the contribution 
of the three-body elastic collisions (cubic in the density) 
to the energy may largely overcome the one arising from 
the two-body terms (quadratic in the density). In this 
situation the condensate appear extremely dilute with 
respect to the two-body collisions but somewhat dense 
with respect to the three-body collisions 14]. 

Theoretical estimates of the three-body coefficient 
were given in 0, 0, 0, 0] . This term can be modeled 
by a quintic nonlincarity in the GPE which for rubidium 
atoms is expected to be attractive |9( . The different types 
of instabilities induced by the quintic term may lead to 
interesting dynamical phenomena in ID BEC similar to 
Bose-Novae effect observed in multidimensional case. On 
the other hand, the presence of these instabilities restricts 
the possibility of BEC manipulations in atomic waveg- 
uides The optical lattice, however, can suppress 
or delay some of these instabilities opening the possibility 
for new types of localized excitations. 

The ID quintic NLS is also used to describe a Bose 
gas with hard-core interactions in the Tonks-Girardeau 
regime. Recent works have shown that this approach de- 
scribes well the ground state prop erties, the collective 
oscillations in a parabolic trap |20j and the dynamics of 
shock waves and dark solitons in the gas [2l|, [23 • On the 
other hand, interference phenomena generated with small 
number of atoms (N = 10) were shown to be not well ac- 
counted by the quintic NLS (the interference patterns 
obtained from this equation are much more pronounced 
than the experimental ones) j2^. For larger number of 
atoms and for weak density modulations, however, the 
description of the nonlinear excitations in the Tonks- 
Girardeau regime by means of a quintic NLS equation 
is expected to be valid. 

The aim of the present paper is to study localized 
states of the quintic nonlinear Schrodinger equation with 
periodic potential for both attractive and repulsive in- 
teractions. To this regard we use a variational method, 
a self-consistent approach and direct numerical integra- 
tions, to show that the presence of the optical lattice al- 
lows to stabilize solitons of the attractive quintic ID NLS 
against collapse or decay. In the case of repulsive inter- 
actions the optical lattice is found to be crucial for the 
existence of stable bright matter waves. The existence 
of localized excitations in this system is shown to be as- 
sociated to the existence of unstable localized solutions 
similar to the Townes soliton of the cubic 2D NLS equa- 
tion (we refer to them as gap-Townes solitons). These 
unstable solutions are found to be located in the forbid- 
den zones of the band structure, very close to band edges 
and separate decaying states from stable gap-solitons. 
The existence curve of these solutions characterizes the 
critical threshold for the occurrence of the delocalizing 
transition and the existence of gap solitons appears to 
be a mechanism for arresting collapse in ID BEC in OL 
with three-body interactions. We also investigate dissi- 
pative effects on gap solitons induced by the imaginary 
part of the three-body interactions. We show that when 



the imaginary part of this interaction is small compared 
to the real part, localized states exist for very long times 
and their behavior is qualitatively similar to the one of 
the undamped case. For moderate and large damping, 
however, all localized states in the band gap decay into 
the Bloch state at the band edge. We find remarkable 
that even in this case a reminiscence of the existence of 
gap-Townes solitons survives in focusing-defocusing cy- 
cles observed when the number of atoms is close to the 
critical value for their existence. We also find that these 
localized states can be stabilized into breather-like excita- 
tions by a linear amplification term modeling the feeding 
of atoms in the condensate from the thermal cloud in the 
presence of nonlinear damping. 

We remark that in absence of periodic potential the 
quintic ID NLS equation has a behavior similar to the 
2D NLS equation with cubic nonlinearity. The interplay 
between dimensionality and nonlinearity has indeed been 
used in the past to investigate collapse in lower dimen- 
sional NLS, the critical condition being D(n — 1) — 4 = 
where n the order of the nonlinearity in the equation and 
D is its dimensionality [24| • From this point of view the 
quintic NLS can be viewed as a ID model for the 2D GP 
with cubic mean field nonlinearity and we expect that 
the results discussed in this paper will apply also to this 
case. 

The paper is organized as follows. In Section II we 
introduce the ID model and discuss the range of appli- 
cability to different physical situations. In Section III we 
study the quintic ID NLS equation with attractive in- 
teractions both in absence and in presence of a periodic 
potential. The effectiveness of the optical lattice to sta- 
bilize localized states of the GPE in presence of two-body 
and three-body interactions is investigated by means of a 
variational analysis and by direct numerical simulations. 
We use a self-consistent method to investigate gap soli- 
tons and band structures. The existence of gap-Townes 
solitons is shown in Section IV and the role of usual gap 
solitons in arresting collapse is discussed. In Section V we 
investigate the case of repulsive interactions by means of 
variational analysis, self-consistent method and numeri- 
cal simulations. The existence of unstable Townes soliton 
with energies above the band edges and their role in the 
delocalizing transition is also investigated. In Section VI 
the dissipative effects introduced by an imaginary part 
of the three-body interaction on gap solitons are stud- 
ied both by a modified variational analysis and by direct 
numerical simulations. The effects of atoms feeding from 
the thermal cloud are also considered. Finally, in the last 
section, the main results of the paper are briefly discussed 
and summarized. 



II. THE MODEL 

Let us consider a BEC with two and three-body in- 
teractions immersed in an optical lattice and an highly 
elongated harmonic trap. In the mean field approxima- 
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tion the system is described by the following 3D GPE 

M 



Equation (J5J is obtained from the following Hamilto- 
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ihu t = - — S/ 2 u+ [—(ujj_p 2 + ujIx 2 ) + V op t{x)]u 

+ gi\u\ 2 u + g 2 \u\ 4 u, (1) 

with p denoting the radial distance, 31,(72 the nonlin- 
ear coefficients corresponding to the two and three-body 
interactions, respectively, lo±,uj x the radial and longitu- 
dinal frequencies of the anisotropic trap (u>± ^S> uj x ) and 
Vopt(x) = Vb sin 2 (kx) the optical lattice applied only in 
the longitudinal direction. The coupling constant of the 
two body interaction is related to the s-wave scattering 
length a s and to the mass m of the atoms by the usual 
relation g\ = Ai:h 2 a s /m. In view of the strong anisotropy 
of the trap we can average the three dimensional inter- 
action over the radial density profile to reduce the prob- 
lem to an effective one dimensional one. More precisely, 
we consider solutions of the form u(r,t) = (j)o(p)'ijj(x,t) 

where <t>o — A / exp(— -^r) is the ground state of the 



radial linear equation 
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Multiplying both sides of the GPE by <fio and integrating 
over the transverse variable, we obtain the following quasi 
1D-GP equation 
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For a discussion on the limits of applicability of this ID 
reduction see |25| . Introducing the dimensionless vari- 
ables 

1 Er u> 2 

t = tv,x — kx. v=——, a=-—r, (4) 

we reduce Eq. 3 to the following normalized ID GPE 
equation with cubic and quintic nonlinearities 
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where a s ~ a s og, with a S Q the constant scattering length, 
g G [—1,1] and Er — h 2 k 2 /2m is the recoil energy 
of the lattice. In this normalization th e wavefun ction 
has been scaled according to ijj — > ijj\j2a S QUj±/v and 
the parameters s,x, are defined as e = Vq/{2Er), x = 
92 {y /2a S QUJi_) 2 jZi^hva^. In the following we will be 
mainly interested in the case 3 = for which it is conve- 
nient to rescale the wavefunction as if) — > ( g^fl Eh ) 1|, V- 
In this case the relation between the dimensionless num- 
ber of atoms N and the physical one N p is N p — 
{^-K 2 a 4 L E R /(g 2 k 2 )) 1 / 2 N (values of N p are typically in the 
range 4000 - 5000 for g 2 /h w 10~ 26 cm 6 s' 1 , a± = 2p m, 
A = 8.23 ■ 10~ 7 m, m = 1.44 • 10~ 25 kg). 



H = j dx[\^ x \ 2 + ^\ 4 + ^\ e + V(xM 2 }, (6) 

as iipt = SH/Sip*. This equation appears also in the con- 
text of nonlinear photonic crystals [26j . The relevance of 
the two-body interactions (cubic term) for the formation 
of localized excitation of soliton type has been largely 
investigated in the past decades in the field of nonlin- 
ear optics and the existence of bistable solitons in optical 
media with Kerr (cubic) and saturable (quintic) nonlin- 
earity was established for the case of a channel waveguide 
with step potential |26| . 

The cubic-quintic NLSE also appears as a model for 
the propagation of optical pulses in double-doped op- 
tical fibers |27| with an effective refraction index given 
by n — n 2 + n^E] 2 , where n 2 = n dl + n d . 2 ,n± = n d2 , 
and n di are the contributions of dopants to the refrac- 
tion index (values and signs can be changed by proper 
choice of dopants). Although Townes solitons have been 
presently found only in connection with collapse of 2D 
elliptic shaped intensive beams in homogeneous focusing 
Kerr media [2^] , we expect gap- Townes solitons to be ob- 
served also in double-doped fibers with Bragg gratings in 
the form of optical pulses and in photonic crystal fibers 
in the form of spatial optical solitons. 

Finally, we remark that a pure quintic NLS equation 
with repulsive interactions also appears in connection 
with a Tonks- Girardeau gas in the local density approx- 
imation [19|. It this context the field equation is 
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also written in normalized form as 
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ax 2 u + ecos(2x)u — \u\ 4 u — 0, (8) 



which corresponds to the case g — 0,\ = 1 in Eq.JSJ). 
Here normalization has been made by introducing di- 
mensionless variables x = kx,t — tu,u — (fry/w/k, with 
J dx\tfi\ 2 — N p = — . In the next sections we use the vari- 
ational approach [? ] , the self consistent method [23] and 
direct numerical integrations, to study localized states of 
the quintic NLS in presence of an optical lattice for both 
attractive and repulsive interactions. 



III. ATTRACTIVE INTERACTIONS 

To obtain analytical predictions for the existence of 
stable localized solutions of Eq.© with a — and for 
attractive interactions {g,x < m Eq.lJ^jl). we use the 
variational approximation (VA) with a Gaussian ansatz 
for the fundamental soliton 
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(9) 



4 




FIG. 1: Early stages of the evolution toward collapse (left 
curves) or decay (right curves) of the critical soliton (central 
curve) in absence of external potentials (a = 0,e = 0). The 
critical soliton in the center refers to the value \ = — 1 an d 
N = 2.7207, while the collapse and the decay evolution refer 
to x — —1-01 and \ = —0.99, respectively, with all other pa- 
rameters unchanged. Snapshots were taken at regular time 
intervals of 5.2 for collapse and of 10 for decay. Plotted quan- 
tities are in normalized units. 



with (i denoting the chemical potential and A, a, the 
amplitude and the square root of the reciprocal width 
of the soliton, respectively. Following the standard VA 
(see [2^ for a review) we derive the following effective 
Lagrangian 
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from which the stationary equations dL /da - 
for parameters a and A are obtained as 
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Here the number of atoms has been expressed in terms 
of a, A, as N = A 2 ^. Note that for s = 0, g = 0, Eqs. 
(|llfl predicts a maximum value of N for the existence of 
soliton given by N cr = \/3 3/2 7r/( 2 x)- This value should 
be compared with the critical norm Ncr — ( 7r /2)\/3/|x| 
obtained for the exact Townes soliton of Eq. J5J with 
e = 0,g = 0,a = 0, 31] 



ip„ = exp(i/it)V>, V 



(3/Vlx|) 1/4 
^/cosh(2 yfjlx) 



(12) 



We see that the VA value deviates from the exact one by 
only 5%. Notice that the above solution has zero energy 
and is marginally stable (dN / d/i = 0) so it exists only for 
the single value of the norm Since for TV > 




I Gap Soliton 
□ .0 0.5 




A 




FIG. 2: Panel a. Lower band and localized bound states 
of the quintic NLS with attractive interactions for parameter 
values x = e = 5. The arrows show the decay pattern 
of the Townes bound state (open circle) for N < N cr (upper 
arrow) and N > Ncr (lower arrow). Panel b. Townes soliton 
(thick line) corresponding to the unstable bound state level 
(open circle) in panel a, with N — N cr — 0.4616 and E = 
—2.1576. Panel c. Gap soliton (thick line and filled circle 
in panel a) obtained from the unstable Townes state for a 
slightly above critical value of N (N = 0.4657, E = -2.2177). 
Panel d. The Bloch state (thick line) at the bottom of the 
band into which the Townes soliton decays for a slightly sub- 
critical value of N (N = 0.4586, E = -2.1531). The OL 
is depicted as a thin dotted line scaled by a factor 100 and 
shifted down by .05 in panel b,c, and scaled by a factor 10 in 
panel d. Plotted quantities are in normalized units. 



Ncf* the solution collapses in a finite time, while for N < 
Nc^) it decays into the uniform background, it appears as 
the analogue of the Townes soliton of the cubic 2D NLS 
equation with attractive interactions. In Fig^we depict 
the early stages of the time evolution toward collapse 
(decay) for an initial norm which is slightly increased 

(decreased) with respect to N$ . 

For e^0 Eq. Ijllll shows that the number of particles 

1 /2 

attains a minimum value -/V t hr — [67rv3(l — 8ee~ 2 )] 
at a = 1/2, i.e. there exists a threshold value of the 
norm which is necessary to create a soliton. The expres- 
sion for iVthr shows that the threshold exists for a rela- 
tively weak lattice and disappears if e exceeds the value 
£0 = e 2 /8 w 0.9236. From this analysis it is clear that 
the optical lattice is quite effective to stabilize solitons 
against decay when N is in the interval N t hi < N < N CI . 
Collapse and decaying into extended states is expected 
for N > N cr and N < N t hr, respectively. The VA makes 
it possible to predict stability of the solitons on the basis 
of the Vakhitov-Kolokolov ( VK) criterion [24| , according 
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FIG. 3: Time evolution of the Townes soliton in Fig|2] (central 
panel) compared with the one obtained from numerical inte- 
gration of the quintic NLS with two slightly different values 
of N: N = l.lJVcr (right panel) and TV = 0.9N cr (left panel). 
Parameters are fixed as in Fig. [5] and plotted quantities are 
in normalized units. 



to which a necessary condition for stability is given by 
dn/dN < 0. Notice that the above variational equations 
are very similar to the ones derived in Ref. [32] for the 
2D cubic NLS (with exactly the same value of eo), this 
being a further confirmation of the analogy between the 
ID quintic NLS and the usual 2D GPE. We remark that, 
since N cr does not dependent on e, the prediction of the 
VA is that the optical lattice does not affect collapse. We 
shall check these results in the next section. 



IV. GAP-TOWNES SOLITONS AND 
ARRESTING COLLAPSE IN ID BEC IN OL 

The existence of the exact Townes soliton l|12|) of the 
quintic NLS with e = makes natural to ask whether 
this solution could exist also in presence of the OL. In 
view of the analogies between ID quintic NLS and 2D 
GPE, the existence of this solution for e ^ could shine 
some light on the role of the OL in controlling collapse, 
as well as, on the origin of the delocalizing transition 
observed in these systems 133]. In the following we use 
a self consistent approach |30| to determine the critical 
value of N for the unstable solution to exist. 

In FigE^ we report the band structure and the local- 
ized states obtained with the self-consistent approach for 
the case g = 0, a = in Eq. J5J. We see that there are 
two bound state levels just below the band. The level 
closer to the band edge (open circle) corresponds to an 




FIG. 4: Left Panel. Delocalizing transition threshold ob- 
tained by direct numerical integrations (filled dots) of the 
quintic NLS with % — —1-0. The horizontal dotted line cor- 
responds to the critical threshold for collapse in absence of 
the OL: iV cr = f v/f- Right Panel. Ener gies versus e of the 

gap soliton (bottom curve with open squares) , Townes soliton 
(middle curve with filled circles) and Bloch state at the bot- 
tom of the band (top curve with open circles) for the quintic 
NLS with x = —1-0. Plotted quantities are in normalized 
unit. 



unstable localized state while the other corresponds to a 
stable gap soliton. The corresponding wavefunctions are 
depicted in Figs. respectively. We find that for a 

fixed value of e there is only one value N for which the 
unstable state exists and its behavior resembles the one 
of the Townes soliton depicted in Fig. ^ In this case, 
however, we have that for a slightly overcritical value of 
N the unstable soliton starts to shrink and its energy de- 
creases as for collapsing solutions until it reaches the gap 
soliton level below, at which the shrinking stops. The 
transition from the unstable Townes soliton to the stable 
gap soliton is represented in Fig|2]} by the lower arrow. 
If the norm is slightly decreased below N cr the unstable 
state completely delocalizes into the Bloch state at the 
bottom of the band shown in Fig[2i (this decay is repre- 
sented in FigEk by the upper arrow). Due to this behav- 
ior we refer to the unstable state as gap- Townes soliton, 
this name being also justified by the fact that for e — > 
the critical norm reduces to N* r (see left panel of Fig. . 
The decaying property of the gap- Townes soliton is illus- 
trated in Fig|3| where the time evolution obtained from 
direct numerical integrations of Eq. |S](with g = 0, a = 0) 
is shown. The left (right) panel of this figure shows the 
decay into the extended state (gap soliton) when the ini- 
tial norm is slightly decreased (increased) with respect to 
N cr . Notice that the transition into the gap soliton state 
is much more rapid than the decay into the extended 
state, this being a reminiscence of the collapse occurring 
at e = 0. Since the energy of a slightly overcritical gap- 
Townes solution should go to — oo, as for any collapsing 
solution, we see that the existence of a stable localized 
state which lies in energy below the unstable state is in- 
deed a mechanism for arresting collapse in the system. 
Also note that a gap- Townes soliton appears as a separa- 



G 



























— j 




k 




n 





20 40 60 80 100 120 
X 



FIG. 5: Time evolution of an intersite symmetric intrinsic 
localized mode for the attractive case x — ~ 1 with N = 2.847 
and e = 0.1. Notice that collapse starts after the mode has 
decayed into an onsite symmetric one (ground state). The 
periodic potential is shown together with the initial state at 
t — 0. Plotted quantities are in normalized unit. 



trix between localized and extended states. The fact that 
its energy lies below the lower band edge implies the exis- 
tence of a delocalizing threshold in the number of atoms 
below which localized state of the quintic GPE with OL 
cannot exist, this being a remarkable difference between 
cubic and quintic nonlinearities (in the cubic case the 
threshold exists only in the higher dimensional case [33|p. 
In Fig. 0Ji the threshold curve separating localized and 
extended states in the (N, e) plane is depicted (this curve 
coincides with the norm of the gap-Townes soliton as a 
function of e). Notice that for e = the norm coincides 
with the exact value ivjr of the solution in 112(1. sug- 
gesting that the family of gap-Townes soliton for e ^ 
originates from the Townes soliton state at e = 0. 

Numerical investigations show that collapse occurs 
when the strength of the optical lattice becomes small 
and the number of atoms overcomes the critical thresh- 
old for Townes soliton at e = 0. For larger values of 
N the width of the soliton becomes small compared to 
the period of the periodic potential and the OL becomes 
ineffective for arresting collapse. We find that solitons ex- 
tended on several potential wells are the ones which are 
better stabilized by the optical lattice against collapse. 
An indication of occurrence of collapse at small values of 
e is obtained from the right panel of Fig. 0] in which we 
depict the energies of the Bloch state at the bottom of the 
band, of the gap-Townes soliton immediately below the 
band and of the gap soliton, are reported as a function of 
e. We see that while the energies of the Bloch state and 
of the gap- Townes soliton monotonically increase with 
decreasing e, the energy of the gap soliton reaches a max- 
imum around e ~ 1 and then starts to decrease as e — » 0. 
This behavior is consistent with the fact that a collapsing 
solutions with energy equal to — oo should exist at e = 0. 
A study of collapse for small e is numerically difficult 



to perform and requires more investigations. It is inter- 
esting to note, however, that besides stable gap-soliton 
the OL allows the existence of metastable states which 
are symmetric around a maximum of the potential (in- 
stead than a minimum as for usual gap-solitons). These 
intersite symmetric states have higher energy than the 
onsite symmetric ones and decay into the ground state 
after some time. Such states can be used to delay the 
time of collapse for N > N CI . Indeed, the collapse of 
these metastable states begins when their instability sets 
in, as one can see from Fig. [S] (notice that the matter 
moves first into a single potential well and then starts to 
collapse) . 



V. REPULSIVE INTERACTIONS 

To investigate the existence of localized states in the 
repulsive quintic NLS with periodic potential we apply 
first a variational analysis based on the following ansatz 
for the soliton waveform Hal 



, , . . sin(aa;) , . 
yj(x, t) = A exp {—ifit) . 



(13) 



Performing the same analysis as before, we get the aver- 
aged Lagrangian 
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from which the following equations for the soliton param- 
eters are obtained 
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Using the fact that the norm and the soliton parameters 
are related by N = it A 2 /a, one can show that 



A = 
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i.e. the VA predicts the existence of a single set of soliton 
parameters a, A for a given norm N and stren gth e of the 
OL. This is similar to what obtained in Ref. |33j for the 
localized solutions of the repulsive 2D GPE with OL. In 
the case of a Tonks- Girardeau gas in OL these parameters 
can be estimated as follows . From Ea. (|l(jfl we get that 
A w (30A^£/227r) 1 / 6 and from the expression of N we 
get for the soliton width uj=l/a» (22A^ 2 /(307r 2 e)) 1 / 3 . 
For a number of particles N = 314 (corresponding in 
physical units to N p = 100) and for Vq = AEr,so e = 2, 
the width of the gap soliton is estimated to be w « 15.5 
corresponding to 3 — 4 lattice periods. 

The band structure associated to gap solitons of the 
repulsive quintic NLS equation has been investigated by 
means of the self-consistent method. In Fig. we show 
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FIG. 6: Panel a. Lower band and localized bound states of the 
quintic NLS with repulsive interactions for parameter values 
X = 1, £ = 2.0. The arrows show the decay pattern of the 
Townes bound state (open circle) for TV < N cr (upper arrow) 
and ./V > N cr (lower arrow). Panel b. Townes soliton (thick 
line) corresponding to the unstable bound state level (open 
circle) in panel a, with TV = TV cr = 1.5565 and E = -0.1047. 
Panel c. Gap soliton (thick line and filled circle in panel a) 
obtained from the unstable Townes state for a slightly above 
critical value of TV (TV = 1.5965, E = 0.2537). Panel d. The 
Bloch state (thick line) at the bottom of the band into which 
the Townes soliton decays for a slightly sub-critical value of 
TV (TV = 1.5065, E = -0.1097). The OL is depicted as a 
thin dotted line scaled by a factor 20 and shifted down by .05 
in panel b,c, and scaled by a factor 10 in panel d. Plotted 
quantities are in normalized units. 



the lower band and the localized bound states for the 
case g = 0, a = in Eq. (J5J. We see that above the 
band edge there are two bound state levels, one imme- 
diately above (open circle) corresponding to an unstable 
gap- Townes soliton, the other more separated from the 
band edge corresponding to a gap soliton. The wavefunc- 
tions of these bound states are depicted in Figs. Eb,c. In 
analogy with the attractive case, we have that the unsta- 
ble soliton exists only for a critical value TV c . r of the norm 
and has a behavior similar to the gap- Townes soliton de- 
scribed before. The decaying property of the repulsive 
gap- Townes soliton is clearly illustrated in Figd where 
we show the time evolution, as obtained from direct nu- 
merical integrations of Eq[S] (with g = 0,a = 0), of the 
gap- Townes soliton in FigEh (central panel). The left 
and right panels of this figure show the decays into the 
extended (Bloch) and localized state (gap soliton) when 
the norm of the initial condition is, respectively, slightly 
decreased or increased with respect to N cr . Thus, also 
in this case a gap-Townes soliton appears as a separa- 
trix between localized and extended states. We remark, 
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FIG. 7: Time evolution of the Townes soliton in FigEJas ob- 
tained from numerical integration of the quintic NLS with re- 
pulsive interaction and three slightly different values of num- 
ber of atoms: TV cr = 1.5565 (central panel), TV = 1.05TV cr 
(right panel) and TV = 0.95TV cr (left panel). Other parame- 
ters are fixed as in Fig|S| Plotted quantities are in normalized 
units. 



however, that in this case the existence of gap-solitons is 
a direct consequence of the band structure (it could not 
exist without the OL). The fact that the energy of the 
Townes soliton is slightly above the upper band edge is 
consistent with the existence of a threshold in the number 
of atoms (given by the critical norm of the gap- Townes 
soliton) below which localized state of the quintic GPE 
with OL cannot exist. In the left panel of Fig. [S] we 
show the delocalizing curve in the plane (TV, e) separat- 
ing localized from extended states. This behavior is very 
similar to what reported for repulsive 2D solitons of the 
GPE in optical lattices |33|. In the right panel of the 
same figure we show the energies of the Bloch states at 
the top of the band, of the Townes soliton immediately 
above the band and of the gap soliton, as a function of 
e. In contrast with the attractive case we see that all 
energies monotonically increase with decreasing e. The 
absence of a maximum in the energy curve of the gap 
soliton at small values of e also indicates the absence of 
collapse in this case (compare the right panels of Figs 0] 
and©. 



VI. GAP-TOWNES SOLITONS IN PRESENCE 
OF DISSIPATION 

In the previous sections we have investigated the gen- 
eral properties of localized states of ID BEC in OL with 
the elastic three-body interactions modeled by a real 
quintic nonlinearity. It is known, however, that the three- 
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FIG. 8: Left Panel. Delocalizing transition threshold ob- 
tained by direct numerical integrations (filled dots) of the 
quintic NLS with \ — 1.0. Right Panel. Energies versus e 
of the gap soliton (top curve with filled triangles), Towncs 
soliton (curve with filled circles) and Bloch state at the top of 
the band (curve with open circles) for the quintic NLS with 
X = 1-0. Plotted quantities are in normalized units. 



body interactions bears also an imaginary component 
corresponding to inelastic collisions which can be mod- 
eled by a damping term of the form —i^\ij}\ A 'il) in the right 
hand side of Eq. (0. Recent theoretical studies [lf^-fl8l] 
estimate for Rb a negative value of the real part of the 
quintic nonlinearity of the order of 10~ 26 — 10~ 27 cm 6 /s 
. From the experiment in Ref. jsij one can deduce the 
imaginary part of the three-body interaction for Rb to 
be of the order 10 _30 cto 6 /s. These numbers imply a 
ratio between the imaginary and the real part of this 
interaction of the order 10~ 3 — 10~ 4 . Other data gives 
values of the damping constant of the order 10 _28 cm 6 /s. 
Due to the uncertainty intrinsic in some of these num- 
bers (presently there are no measurements of the real 
part of the three-body interactions) in the following we 
take the damping constant 7 to be a free parameter and 
investigate both the underdamped and the overdamped 
regimes. 

The instability of gap-Townes solitons against small 
fluctuations in the number of the atoms rises the ques- 
tion of how the scenario of the previous sections changes 
in presence of dissipation. To answer this question we 
concentrate on the case of attractive quintic interactions 
(qualitatively similar results holds also for repulsive inter- 
actions) in presence of a quintic dissipation. The results 
are displayed in Fig. [5| In the top three panels of this 
figure we depict the time evolution of the Townes soliton 
in Fig. in presence of a small damping. We see that 
in spite of a slow broadening of the Townes soliton (mid- 
dle panel), the situation remains qualitatively similar to 
the undamped case depicted in Fig. [3] Notice that for 
an overcritical number of atoms (right top panel) sev- 
eral focusing-defocusing cycles between the gap-Townes 
soliton and the gap soliton occur (these oscillations are 
induced by the slightly overcritical starting value of N). 
From the left panel of this figure we also see that for an 
undercritical number of atoms the decay into the Bloch 
states occurs faster than for the zero damping case. At 



very long times the gap-Townes soliton in the middle 
panel of Fig. eventually decays into the Bloch stated 
at the edge of the band. 

For larger damping the situation is different, as one 
can see from the bottom three panels of FigO Notice 
that for an overcritical value of N the gap-Townes soliton 
(right panel) decays into an extended state after perform- 
ing only one focusing-defocusing cycle. In this case the 
dissipation completely suppresses collapse and localized 
states becomes unstable against decay into the lowest en- 
ergy Bloch state. It is remarkable, however, that a signa- 
ture of the existence of gap-Townes soliton remains even 
for moderately large dampings in the focusing defocusing 
cycles observed when the decreasing (due to dissipation) 
number of atoms crosses the critical value for existence of 
gap-Townes soliton. This fact could indeed be used in ex- 
periments to detect the existence of gap-Townes solitons 
in presence of damping. 

An analytical explanation of the existence of focusing- 
defocusing cycles can be obtained by means of a mod- 
ified variational analysis [sM l8rl | in which the damping 
is treated as a small perturbation. In this approach we 
assume a time dependent localized state of the form 



4>{x, t) = Aexp(— - — - + ibx 2 + i 
2w z 



(17) 



with parameters 77; = {A, w, &,</>} also time dependent. 
The equations of motion for these variables are obtained 
from the averaged unperturbed Lagrangian of Eq.JSJ 
with 3 = 



1 



L = -^A 2 w[w 2 bt+2(l) t +A;+4w 2 b 2 -2ee- w 



2X 
3\/3 



A 



as 



d dL dL 



dt drji 



dL f 

-- — = dx{R4,* m + R*^ m ) (18) 
t or]i 7 -00 



where R = — ij\ip\ 4 ip, and the subscript T)i denotes the 
derivative with respect to soliton parameters. From the 
above generalized Euler-Lagrange equation one derives 
the following equations for the number of atoms N and 
soliton width w 



N 3 

N t = -3 7 /3— 

wtt = 4r(l + xPN 2 ) - 8ewe~ 



(19) 



with (3 = 2/(3 v / 3¥). 

In Fig ^]| we show the time evolution of the soliton 
width for different values of the dissipation constant. 
We wee that for very low dampings the width oscillates 
around a constant value, meaning that the localized state 
undergoes periodic focusing-defocusing cycles in agree- 
ment with the right top panel of Fig. By increasing 
the damping we see that after a certain number of cycles, 
the soliton width starts to grow monotonically signaling 
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FIG. 10: Left panel. Time evolution of the width of the 
Townes soliton as obtained from Eqs. I19l with e = 2, \ = — 1 
and damping parameters increased by a factor 10 going from 
the bottom curve (7 = 0.0001) to the top one (7 = 0.1). 
The last three curves have been shifted respectively by 2, 4, 6, 
to avoid overlapping. The initial value of atoms was detuned 
from the critical value N = 1.254 corresponding to the unper- 
turbed Townes by 0.05. Right panel. Time evolution of the 
soliton width obtained from the damped quintic GPE (contin- 
uous line) and from the variational eauations ll9l (dotted line). 
Parameter values are taken as e = 2, x = — 1> 7 — 0.01. The 
comparison with the variational results was done by rescaling 
the soliton width obtained from the GPE by a factor 11.5 
to have the same initial condition. Plotted quantities are in 
normalized units. 



FIG. 9: Time evolution of Townes soliton in Fig|5|as obtained 
from numerical integration of the dissipative quintic NLS with 
attractive interaction. The top three panels refer to three 
slightly different values of number of atoms: N cr = 0.87248 
(central panel), N = l.liV cr (right panel) and N = 0.9iV cr 
(left panel) with a dissipation parameter 7 = 0.00025. The 
bottom three panels are the same but for 7 = 0.01. Other 
parameters are fixed asx = — l,e = 5. Plotted quantities are 
in normalized units. 
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the decay of the gap soliton into an extended state. No- 
tice that even for relatively high dampings there is one 
focusing-defocusing cycle, this being in good qualitative 
agreement with what observed in the bottom right panel 
of Fig. [5] In the right panel of Fig. ^|we make a quan- 
titative comparison between the dynamics of the soliton 
width as obtained from the VA and from the damped 
quintic NLS. From this figure we see that although the 
main feature of the phenomenon is correctly captured 
(i.e. existence of only one focusing-defocusing cycle be- 
fore decay), the variational equations are not in good 
quantitative agreement with the numerical integration 
of the full system. This may be due to the fact that 
the ansatz in Ea. i|17fl is more appropriate for strongly 
localized states than for extended gap- Townes ones. It 
is interesting, however, that for e = the variational 
equations become similar in form to the ones considered 
for the 2D cubic NLSE in presence of a cubic dissipa- 
tion — ij\ip\ 2 ip In this case one can show that the 
equation for w reduces to a modified Airy equation for 
which a focusing-defocusing cycle always exists. When 



FIG. 11: Left panel. Time evolution of the gap-Townes soli- 
ton in the lower right panel of Fig. |5] in presence of quin- 
tic dissipation 7 = 0.01 and linear amplification of strength 
8 = 8.95 ■ 10~ 4 . Right panel. Time evolution of the gap- 
Townes soliton amplitude in the left panel. Plotted quantities 
are in normalized units. 



e 7^ 0, Eq.s l|19|) predict that for decreasing 7 the number 
of focusing-defocusing cycles grows, this being in good 
agreement with direct numerical simulations. Moreover, 
from the left panels of Figs. 21 an d El we see that the 
critical value of the number of atoms for existence of 
gap-Townes soliton decreases by increasing e, this mean- 
ing that in strong optical lattices the effect of the quintic 
damping on gap-Townes solitons is effectively reduced 
(due to the reduced norm), this being confirmed also by 
numerical simulations. 

We remark that besides the decay phenomena induced 
by the quintic dissipation one could also consider the 
feeding of atoms in the condensate from the thermal 
cloud 38]. The injection of matter in the condensate can 
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be modeled with a linear imaginary term (of sign opposite 
to 7) in the NLS equation of the form iSip. This term can 
balance the damping and stabilize localized states against 
decay even in the overdamped regime, this leading to the 
formation of dissipative gap solitons. We illustrate this 
phenomenon in Fig. 1111 where the time evolution of the 
slightly overcritical gap- Town soliton shown in the bot- 
tom right panel of Fig. [!|| is reported for a linear am- 
plification term of strength 6 = 8.95 • 10~ 4 . From this 
figure we see that the localized state instead of decaying 
into the Bloch state at the bottom of the band (see Fig. 
|5J), it executes breather-like oscillations (i.e. focusing- 
defocusing cycles) with emission of radiation. In terms 
of band picture this breathing state can be thought as a 
non stationary states with energy oscillating between the 
gap-Townes soliton and the gap soliton in Fig. [2^. We 
believe this being a general mechanism for existence of 
breather-like excitations in ID BEC in OL in presence of 
damping and amplification. 

VII. CONCLUSIONS 

In this paper we have used the periodic ID NLS equa- 
tion with quintic nonlinearity to investigate localized 
states in ID BEC in OL with three-body interactions. 
The existence of unstable solitons similar to the usual 
Townes soliton of the cubic 2D NLS has been established. 
These states have been shown to have energy located in 
the forbidden zones of the band structure, very close to 
band edges, separating decaying states from stable lo- 
calized ones (gap-solitons). The existence of gap soli- 
tons appears to be a mechanism for arresting collapse in 
attractive low dimensional BEC with three-body inter- 
actions in OL and the analysis of gap-Townes solitons 
appears to be important for characterizing the delocaliz- 
ing transition in these systems. To this regard we remark 
that the region of existence of localized states in the N, e 
plane is bounded by two limiting curves. The first curve 
separates extended states from localized ones and there- 
fore characterizes the delocalizing transition. This curve 
coincides with the curve of existence of gap-Towenes soli- 
tons. For the ID quintic NLS with periodic potential 
the delocalizing transition exists for both attractive and 
repulsive interactions. Delocalizing curves have been in- 
vestigated in Fig.s 01 |H1 A similar behavior is expected 



to be valid also for the standard CPE with periodic po- 
tential in higher dimensions. The second curve separates 
the localized solutions from the collapsing ones, giving 
an upper bound for the existence of gap-solitons. This 
curve obviously exists only for attractive interactions (for 
repulsive interactions a negative effective mass, although 
producing an attraction among atoms which leads to the 
formation of gap solitons, is usually insufficient to induce 
collapse). Due to the numerical difficulties involved in 
the study of collapsing solution, the threshold curves for 
collapse have not been investigated in this paper. For 
this it will be useful to develop an analytical approach to 
provide estimates of the critical values for collapse [3!| ■ 

We have also investigated the influence of a dissipa- 
tive imaginary part of the three-body interaction on gap- 
Townes solitons. In particular, we have shown that the 
damping introduces effects which are important when 
the ratio between the imaginary and the real part of the 
three-body interaction is not small. For the case of ru- 
bidium this ratio is expected to be small and the situa- 
tion described for the underdamped case should be qual- 
itatively valid. It is remarkable, however, that even for 
moderate large damping a signature of the existence of 
gap-Townes solitons remains in the focusing-defocusing 
cycles observed when the number of atoms crosses the 
critical value for their existence. 

Finally, we have shown that the presence of a small 
linear amplification term into the damped quintic NLS 
equation, modelling the feeding of atoms from the ther- 
mal cloud, can stabilize gap solitons into breather-like 
excitations even in presence of large damping. 

The results of this paper suggest the possibility to ex- 
perimentally observe gap-Townes solitons in low dimen- 
sional BEC. 
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